Concentration and purification of entanglement for qubit systems with ancillary 

cavity fields 
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We propose schemes for entanglement concentration and purification for qubit systems encoded in 
flying atomic pairs. We use a cavity-quantum electrodynamics setting as the paradigmatic scenario 
within which our proposals can be implemented. Maximally entangled pure states of qubits can be 
produced as a result of our protocols. In particular, the concentration protocol yields Bell states 
with the largest achievable theoretical probability while the purification scheme produces arbitrarily 
pure Bell states. The requirements for the implementation of these protocols are modest, within the 
state of the art, and we address all necessary steps in two specific set-ups based on experimentally 
mature microwave technology. 

PACS numbers: 



I. INTRODUCTION 

Besides the attention that modern physics pays to en- 
tanglement in virtue of the role it plays in testing funda- 
mental properties of quantum systems, there is an enor- 
mous interest in developing strategies and protocols al- 
lowing for the manipulation of these intimately quantum 
mechanical correlations between discrete as well as con- 
tinuous variable systems. Even though its role in compu- 
tational processes is not yet fully understood and is the 
subject of intriguing investigations, it is nowadays rec- 
ognized that the ability of manipulating entanglement is 
crucial for many tasks in quantum communication and, 
in general, quantum information processing (QIP). 

Over the past two decades, we have witnessed many 
proposals for the generation of entangled pairs of quan- 
tum systems, in various settings For instance, it has 
become a routine procedure in linear optics QIP with 
discrete variables, in which parametric down conversion 
is used @ and in its continuous variable version, where 
two-mode squeezed states are effortlessly produced by 
non-degenerate optical parametric amplifiers. In differ- 
ent contexts, entanglement of pairs of atoms [H, trapped 
ions Q, superconducting devices [BJ and atom-photon 
systems @ has already been demonstrated. 

Regardless of the efficiency with which it can be gen- 
erated, entanglement is undoubtedly a fragile resource, 
susceptible to the detrimental effects of communication 
channels and environmental interactions that result in 
a degradation of the desired "quality" of the quantum 
resource. This is particularly true in a scenario of dis- 
tributed QIP, where the exchange of high-quality entan- 
gled states (which should be as pure as possible and car- 
rying a full ebit of quantum correlations) is necessary 
in order to realize perfectly efficient computational pro- 
tocols. This has lead to the development of strategies 
to increase the entanglement and purity of mixed, non- 
maximally entangled states @, [!]■ So far, experimental 
realizations of entanglement purification schemes have 
been limited to linear optics [§, [T3, EH and trapped- 
ions (T^ |. Though lacking such a demonstration, mi- 



crowave cavity-quantum electrodynamics (cQED) holds 
great promise as a medium for QIP. Strong coupling be- 
tween atoms and field modes allows for the creation [l3| 
and distribution [141 ] of entanglement, and the perfor- 
mance of quantum logic operations [l5| |. The new high 
finesse cavities [l6j demonstrated by Kuhr et al. im- 
prove the feasibility a system of distributed QIP employ- 
ing atoms as flying qubits, and cavity field modes as rela- 
tively long-lived static qubits. This system would clearly 
be enhanced by methods for improving degraded entan- 
glement. 

Suggestions have been made for cQED purification 
schemes that recover entanglement only from decoher- 
ence due to cavity decay [l7j, or that rely on repeated 
interactions of the cavity modes with a single flying 
atom Here, we propose schemes for the concentra- 
tion of pure-state entanglement, and the distillation of 
pure-state entanglement from Werner states (assuming 
a sufficiently large fraction of any Bell state) and max- 
imally entangled mixed states (MEMS) [lj|. Our pro- 
tocols rely only on single-atom operations and measure- 
ments, and could be performed using current technology. 

The paper is organized as follows. In Sec. [Til wc in- 
troduce the interaction model used throughout this in- 
vestigation and describe the protocol for the concentra- 
tion of entanglement that represents the first part of 
our study. This probabilistic entanglement concentra- 
tion is accompanied by a simultaneous interconversion 
from flying-qubit to stationary-qubit entanglement. We 
prove that a sufficiently long interaction time allows for 
the concentration of entanglement to the maximal value 
of a full ebit with the highest probability being theoreti- 
cally allowed. The scheme requires resources that can be 
prepared off-line and is extremely close to ideal efficiency 
for values of the experimentally relevant parameters at- 
tainable with present-day technology. In Sec. IIIII we an- 
alyze an entanglement purification protocol that, when 
operating on mixed Bell-diagonal states, extracts a pro- 
gressively more entangled and pure two-qubit state. The 
scheme is probabilistic and requires information gained 
on the population of the cavity fields. We demonstrate 
the wide applicability of our scheme by addressing the 
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purification of two interesting classes of mixed states, 
namely Werner states and maximally entangled mixed 
states of two qubits. Sec. IIVI is devoted to the analysis 
of the feasibility of the schemes we propose in microwave 
cavity- and circuit-QED set-ups. We highlight how the 
recent achievements in this field allow for a foreseeable 
implementation of our proposals. Moreover, we discuss 
the possibilities offered by superconducting devices in- 
tegrated in planar resonators as alternative scenarios in 
which our protocols can be realized. Finally, Sec.fVlsum- 
marizes our results. 



II. CONCENTRATION SCHEME 

Many problems involving the interaction of spin-like 
two-level systems with single-mode bosonic systems can 
be modelled with an effective dipole-coupling formalism 
in which a spin operator (generally proportional to the a x 
Pauli spin operator) is coupled to the electric (or the mag- 
netic) part of an electromagnetic field. This is, evidently, 
the case for neutral atoms or quantum dots coupled to 
optical or microwave fields [1, H3, Hl|. On the other 
hand, this description holds also for a system consisting 
of a Cooper-pair box (in a superconducting-quantum- 
interference device (SQUID) configuration and in the 
charge regime (22| ) integrated into a planar stripline res- 
onator 23], a setting we generally refer to as circuit-QED. 
At the charge degeneracy point, an effective dipole mo- 
ment operator for the SQUID can be written, whose am- 
plitude is proportional to the excess charge in the SQUID 
island [23]. In what follows, in order to fix the ideas 
and introduce the general formalism employed through- 
out our study, we shall use language and terminology 
that are typical of cavity-QED, and we refer explicitly to 
a scheme of Rydberg atoms interacting with microwave 
cavities. In Sec. |IV] we extensively assess the details of 
oimplementation in both cavity- and circuit-QED. 

Let us consider a two-level atom with ground and ex- 
cited states \g) and |e), respectively. The corresponding 
transition frequency is labelled oj . The atom interacts 
with the field mode of a cavity of frequency o>/, described 
by the bosonic annihilation (creation) operator a (a 1 *). 
Within the dipole-coupling interaction assumed here, the 
total Hamiltonian of this atom-field system is 

H = hw a z + hajfa^a + HA (&+ + <x_) (a + a^) , (1) 

where <r+ = |e) {g\ and o~ — \g) (e\ are the atomic ladder 
operators and A is the coupling strength of the interac- 
tion. Employing the rotating wave approximation [24j, 
assuming atom-field resonance (i.e. ljq ~ w/) and enter- 
ing the interaction picture with respect to the free energy 
Hq = hw a z + fno a) a, the coupling reduces to the stan- 
dard Jaynes-Cummings model [251 ] 

Hi nt = h\(a+a + (t-a 1 ) . (2) 



The dynamics arising from this Hamiltonian |J5J| are de- 
scribed by the following time-evolution operator, written 
in the atomic basis {|e), \g)} [26| 

cos (Xty/ n + 1) — zd sin ^^f™^ 
— la 1 — , . , — '- cos(Atvn) 

where n = is the photon- number operator of the 
cavity field. 

Even though this simple interaction model has been 
extensively studied in recent years, the large variety of 
effects for which it is responsible and its vast applicabil- 
ity to problems of QIP make it worth our consideration. 
Here, we give a simple, yet interesting example of its ap- 
plicability to the problem of entanglement concentration. 

Consider a non-maximally entangled pure state of two 
atoms, labelled 1 and 2: 

\^o)u=ot\ee) 12 +/3\gg)i2- (4) 

where a is real and /3 = y/l — a 2 . The assumption of 
this form of the input state does not affect the general- 
ity of our approach. (Other input states may be used if 
converted into the form ([?]) , which can be achieved by 
using Ramsey zones to perform single-qubit rotations.) 
This state can be created using a variant of the scheme 
suggested in Rcfs. 27] , which is able to establish an ar- 
bitrary amount of entanglement in the pure state of two 
preparatory cavities cl and c2, i.e. a state of the form 
cos(7rr)|01) c i. C 2-|-sin(7rT)|10) c i. C 2. Atoms 1 and 2 are pre- 
pared in their ground states, and each is passed through 
one of the two cavities. By arranging a quarter-Rabi cy- 
cle between atom j and the field of cavity cj (j = 1,2), 
and using a single Ramsey zone [l3j to perform a a z rota- 
tion of angle it on one of the atoms, we will transfer the 
entanglement from the cavities to the atoms. Alterna- 
tively, one can use a single cavity sequentially crossed by 
atoms 1 and 2 as in Ref. [2l| . Of course, the preparation 
of the input state can be considered as an off-line step 
that should not be accounted among the requirements of 
the scheme we address. 

We present two versions of our concentration protocol; 
a symmetric scheme, in which a can take any value in 
the range [—1,1], and an asymmetric scheme in which 
we require a > /3. 

A. Symmetric Concentration Scheme 

We allow qubits 1 and 2 to interact with the field of 
cavities a and b respectively, which are spatially sepa- 
rated with each prepared in a single-photon state. A 
sketch of the described setup is presented in Fig. [T] 
The initial state of the overall system is, therefore, 
|*o)i2 <8> |H)a6 which evolves according to the global 
time-propagator i7i a (ii) £g> L^bfe), where Uj^(tj) is given 



by the matrix expression Eq. ([3]) with the replacement 
a, — > cij, at (J = a, 6) and atomic label /x = 1,2. 

For the sake of simplicity, we set the interaction time 
t\ = t<i = t, which results in the evolved state 

|*a/)l2a& = ^[|ee, lp ee )l2ab ~ \99, 1pgg)l2ab ,^ 
- l\eg, ^eg)l2ab ~ Age, ^ge)l2ab}- 

In this equation, we have introduced the normalized two- 
mode states 

\i>ee)ab = N ee \oc COS 2 (V2Ai)|ll)afc-/?sin 2 (At)|00) o6 ], 

\fpeg)ab = A/" eff [a COS (2V2At)|12) afc + /3sin (2At)|01) a6 ], 

\lpge)ab = A/" ge [a COS (2 V2At) |21) o6 + p sin (2At) 1 10) o6 ] , 

|V>g 9 >«f> - AA 99 [asin 2 (V2At)|22) ab - /3cos 2 (Ai)|ll) a6 ] 

(6) 

For the sake of mathematical convenience, we shall con- 
sider the case of < a < (3. These equations show that 
the information initially contained in the input two-atom 
state has been transferred to cavity fields in the sub- 
space {|0),-, |2)j} of the relevant infinite-dimensional 
Hilbert space. To complete the entanglement concentra- 
tion process, we measure the state of each atom in the 
cr z -eigenbasis {|e),|g)}. After postselecting the events 
corresponding to both atoms being found in |e), this op- 
eration results in the projection of the cavity fields onto 
\?Pee}ab with probability 

P(\ip ee )) = a 2 [cos(%/2A£)] 4 + /3 2 [sin (At)] 4 . (7) 

If the interaction time was selected such that 

a[cos(V2At)] 2 = /3[sin(Ai)] 2 (8) 

we produce a two-mode state with a full ebit of en- 
tanglement. Therefore, we have probabilistically con- 
centrated the initial flying-qubit entanglement by trans- 
ferring it to two stationary qubits, here embodied by 
the cavity field modes. This full ebit can, with unit 
probability, be transferred to two flying atoms by ar- 
ranging local interactions with the cavity fields for a 
time such that the transformations \g, 0) a t om ,cavity —> 

\di ^} atom, cavity ) \g, atom, cavity * \e, 0} atom, cavity are re- 
alized. 

Let us now consider whether or not this concentration 
scheme is ideal. To this end, we use the tool provided 
by the entanglement of single pair purification (ESPP), 
Es [28( | . This quantity, which is shown to be a measure 
of entanglement in Ref . 28] , is defined as the maximum 
probability with which a non-maximally entangled qubit 
pair can be converted to a full-ebit state using local oper- 
ations and classical communication (LOCC). It is equal 
to 2|c s | 2 , where c s is the smallest Schmidt coefficient in 
the state under investigation [28[. The input state (U]) 
has Es = 2a 2 . The optimality of the concentration pro- 
tocol is evaluated as follows: As ESPP is a good measure 
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FIG. 1: Entanglement concentration scheme. Two cavities 
act as stationary qubits interacting with flying atomic pairs 
previously prepared in a non-maximally entangled state. Lo- 
cal cavity-atom interactions followed by atomic-state detec- 
tion probabilisticly transfer one ebit of entanglement to the 
cavity field modes. For properly chosen interaction times, 
the probability of obtaining an ebit is arbitrarily close to the 
optimal theoretical value. 



of entanglement, it cannot increase, on average, under 
application of LOCC (29|. Therefore, the protocol for 
entanglement concentration is said to be optimal if it 
conserves ESPP. 

As before the measurement of the atomic pair we have 
an a priori ensemble of states given by Eqs. ([5]), each 
obtained with a precise probability, the object of our in- 
terest will be the average entanglement of the ensemble. 
After measuring the state of the atomic pair, the cavity 
modes have the average ESPP 

(E s ) f = P (|Vw» E S {\i> ee ))+p{\^ eg ))E s (\^ eg )) 

+ p(\^ ae ))E S {\^g e ))+P^gg))E S ^ gg )), 

where p (IV'ij)) is the probability of projecting the field 
modes onto the state = e,g), which has ESPP 

Let us momentarily impose the following condition on 
the interaction time 

cos 2 (XtV2) = 1 — > Xt = kn/V2 (k € Z). (10) 

It is clear that this will set Es (\ip eg )) = Es Qi/ige)) — 
Es (IV'gg)) = 0- We now make a simple trigonometric 
observation: As kn/y2 and kn are incommensurable, 
by letting k grow, the value of sin(fc7r/v2) will fluctuate 
within the range [—1,1], such that the entire range will 
eventually be covered. This is shown pictorially by the 
trigonometric circle in Fig. [5] (a). 

Thus, having set the constraint in Eq. (fTT)|) , we can still 
find a value for the rescaled interaction time Xt such that 
any additional condition we choose to set on sin(Ai) will 
also be satisfied. The only requirement is having a suffi- 
ciently large range of integers k, i.e., a large upper bound 
on the interaction time t [30l ]. For instance, if we select 
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FIG. 2: (a): Distribution on the trigonometric circle of the 
values of sin 8 and cos 8 associated with 8 = kn/\/2 for 
k = 0, ..,50. (b): Distribution of sin 6 1 against k compared 
with the value of = 0.512418, arbitrarily taken from a 

Gaussian ensemble of 100 values chosen in the range [0, l/y/2]. 

s'm(Xt) = 0, we see that the average ESPP of the cavity 
modes is zero. If, on the other hand, sin 2 (Ai) = a/ [3, we 
find that condition © is satisfied, so that Es (\ifj ee )) = 1, 
and (Es)j — P(lV'ee)) = 2a 2 . Thus, for general values 
oft, (E s ) f < 2a 2 , the upper bound being attained when 
both conditions ([5]) and p0[) are satisfied. In this case, 
the procedure gives the maximum probability of produc- 
ing a full-ebit state, and can be said to be optimal. Note 
that in this case, a single round of the concentration pro- 
cedure is all that is necessary to (probabilistically) pro- 
duce a maximally entangled state from any input pair. 

As an example, in Fig. [2] we s how the case of a = 
0.253964, corresponding to y^/ft = 0.512418 and k free 
to span the first 50 integers. Panel (a) clearly shows 
how these fifty values are uniformly distributed over the 
trigonometric circle. On the other hand, in panel (b) 
we compare the distribution of | sin(/c7r/v / 2)| with the 
chosen value for \J a/ ft, (which has been extracted from 
a Gaussian ensemble of values of a lying in the range 
[0,1/^2]). The case is rather favorable as, for k = 4, 
sin(2v / 27r) = 0.513288 « y/a/P- Clearly, for other val- 
ues of a, larger values of k may be required in order for 
conditions ([5]) and (fTT)|) to simultaneously hold. 

However, in practical realizations, the interaction time 
is limited by the coherence time of the system, which is 
given by 

Tcohe 

= min(T at0 m,T field), where 
is the cavity field (excited atomic state) lifetime. Our 
unitary treatment for entanglement concentration has to 
be consistent with the requirement t < r co / le . We there- 
fore change our approach by firstly requiring that con- 
dition be fulfilled, so that all successful runs pro- 
duce a Bell state output. Secondly, the probability of 
success is maximized within some upper bound on Xt. 
Fig. [3] shows this maximum probability, P ma x , against 
the ESPP of the input state Eq. 0$ when the rescaled 
interaction time Xt is limited by 7r, 57r and 207T (dashed, 
full and dotted line respectively). The implications of 
these upper bounds in terms of requirements for r co h e 
are discussed in Sec. IIV1 The analysis of Fig. [3] confirms 
that a sufficiently large upper bound to the interaction 
time brings the protocol to optimality. The maximum 
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FIG. 3: Maximum probability, Pmax, of the entanglement 
concentration protocol yielding a Bell state, plotted against 
the input ESPP (Es) for the case of Xt bounded by tt, 5tt 
and 20-7T (dashed, solid and dotted line respectively). The 
larger the upper bound to the rescaled interaction time, the 
closer the protocol to the ideal performance Pmax = Es, im- 
plying conservation of the initial ESPP under application of 
the concentration protocol [281 ] . 



probability of concentrating entanglement into a full-ebit 
state becomes equal to 2a 2 , and so the ESPP of the in- 
put state is, therefore conserved. Smaller upper bounds 
will, in general, reduce P ma x- Nevertheless, it is impor- 
tant to stress that even at very short interaction times 
(i.e Xt < tt, dashed line), there are input states | v I / o)i2 
allowing for optimal concentration. 

The successful implementation of this concentration 
protocol clearly relies on our previous knowledge of the 
amplitude a in the initial entangled atomic state 1^0)12 
in Eq. (g]). Here, we want to assess the effect of an un- 
certainty in the value of a on probability of success of 
the protocol, i.e., the probability of projecting both an- 
cillary atoms onto their excited state. Our assumption 
is that a follows a Gaussian distribution centered at a 
with standard deviation a . This is perfectly reasonable 
as it is consistent with the protocol we suggested in or- 
der to generate input state (g]), where an uncertainty in 
the velocity of the preparatory atom (usually assumed to 
follow a Gaussian distribution) may result in some fluc- 
tuations in the degree of entanglement between atoms 
1 and 2. The quantity of interest is the mean success 
probability P(\ip) ee ), as given in Eq. (JT]), averaged over 

the distribution G(a,a) = . 1 expf— ^ a ~^ 1. A nu- 
merical analysis reveals that the values of the rescaled 
interaction time giving the maximum success probability 
are unaffected by a, but the success probability decreases 
as the uncertainty in a increases, as should be expected. 
Clearly, the average entanglement of output states will 
decrease with increasing a. 



B. Asymmetric Concentration Scheme 

This protocol is similar to that described above, the 
difference being that the cavities are prepared in the vac- 
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uum state, which eases the implementation of the proto- 
col. After the atoms, prepared in state ([!]), interact with 
the cavities for a time t, the state of the system reads 

\psi{t)) abl2 = [-asin 2 (Ai)|ll) ab + /3|00) afc ]| ff5 )i 2 

-~ sin(2At)[|01) afc |e< ? )i 2 + |10) a6 | 5 e)i 2 ] 
+acos 2 (Ai)|00> Q6 |ee)i 2 (11) 

As opposed to the symmetric scheme, we here require 
a > f3 and set sin 2 (At) = (3 /a. By projecting both atoms 
onto their ground states, we obtain, with the maximum 
allowed probability 2/3 2 , a maximally entangled state of 
the cavity fields. Note that we have imposed only one 
condition on t] thus the optimal interaction time will 
always exist within the limit Xt < ir/2. The required 
form of the input state can be generated from any input 
by means of local rotations performed on qubits 1 and 2 
at properly set Ramsey zones 0, [H[ . 




FIG. 4: Entanglement purification scheme. Two pairs of fly- 
ing qubits (1&2 and 3&4) are prepared in the same mixed 
(Bell-diagonal) state and interact with two cavities (a and b 
respectively). The atomic-state measurement of pair 1&2 ac- 
companied by the projection of the cavity fields onto photon- 
number states realizes a conditional purification of the entan- 
gled mixed state shared by qubits 3&4. 



III. PURIFICATION SCHEME 

We now move on to the description of a simple quan- 
tum state purification scheme 0, HI that is able to oper- 
ate on any Bell-diagonal mixed state. A general remark 
about the tools used in this Section is due. As we mainly 
deal with bipartite mixed states of two qubits, we shall 
quantify quantum correlations in terms of the negativity 
Eji/ [31| . This quantity is strictly related to the positive 
partial transposition (PPT) criterion for the separability 
of quantum states [32]. The partially transposed den- 
sity matrix g is obtained from any given bipartite quan- 
tum state by transposing the variables of only one of the 
two subsystems. The PPT criterion then simply reads 
g > 0, and is necessary and sufficient for the separability 
of any bipartite two-qubit quantum state. The 'nega- 
tivity' E^j- [3l[ , which can be easily determined from the 
density matrix, is defined as the absolute value of the sum 
of the negative eigenvalues of g and directly quantifies the 
violation of the PPT criterion. It is monotonically related 
to the logarithmic negativity, which operatively quanti- 
fies an upper bound to the distillable entanglement. 

Let us consider a state of the form 

j=± 

where |$±) = (l/v^)(|00) ± |11» and |tf±) = 
(1/V2)(|01) ± 1 10)) are the elements of the Bell basis 
and Ylj=±i-Aj + B,j) = 1. Considering that the par- 
tial transposition of a two-qubit density matrix has only 
one negative eigenvalue [l9j , it is straightforward to show 
that if A+ > 1/2, E^ = 2A + - 1. Obviously, Ej^ = 
for an equal mixture of Bell states. 

Let us begin our analysis by assuming the special case 
of A_ = _B_ = so that, imposing normalization, A+ = 



1 — B+ = V . The resulting state reads 

= + (1 - V) |* + )<* + |, (13) 

for which E$ = 2V - 1. The task of our study is the 
purification of p Q by using two replicas of the state and 
an ancillary system. In our scheme, the replicas of p Q 
are encoded in the state of qubits 1 & 2 and 3 & 4, while 
the ancillary system is embodied by two cavities (labeled 
a and b) analogous to those considered in the previous 
Section. In order to illustrate the main features of our 
protocol, we momentarily assume that the cavities are 
prepared in an n-photon state, with n 3> 1. We will 
relax this assumption later on. The initial state of the 
system is thus p ,i2 <8> p ,34 ® \nn) a b{nn\. 

As shown in Fig. 01 one atom from each pair passes 
through each cavity (qubits 1 & 3 pass through a while 
qubits 2 & 4 pass through b) and the interaction time t 
is selected so that y/nXt = tt/2. In order to understand 
the main features of the scheme we propose, it is enough 
to consider just the dynamics within cavity a, where, by 
neglecting the small difference between the effective Rabi 
frequencies y/nXt and y/n+lXt, the following transfor- 
mations occur: 

|ee>i 3 |n> -> |5ff)i3|" + 2> a , \eg) 13 \n) a -> \ge) 13 \n) a , 

\gg)i3\ n )a -> | ee) 13 1 ra - 2) a , |ge)i 3 |n) -> |eg)i 3 |n) a . 

(14) 

The cavity fields are now projected onto \nn) a b with a 
success probability [P 2 + (l — "P) 2 ]/2. Only those terms in 
which both atoms passing through each cavity are in dif- 
ferent states survive the projective measurements. Thus, 
if both atom pairs are initially in the same Bell state, 
the output will be a GHZ-like state of four atoms, while 
cross terms will be annihilated; 
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|<J> + )l 2 |$+) 3 4 
|*+)l2|* + ) 3 4 0, |* + )l 2 |*+)34 



1234 = -j={\ggee) + \eegg)) 123i , 



|*+)i 2 |*+) 3 4 |* 2 )i234 = -^(\9eeg) + \egge)) 123 4. 

(15) 

The effect of the ancillary systems embodied by the 
cavity fields is to filter the right form of correlations 
among the qubit pairs, leaving us with the mixed state 

(1) _ ^ 2 |$ 2 )l234($ 2 | + (1 - Vf |fr 2 )l234(fr 2 | ri _. 
Pl234 ~ V 2 + {1 _ V f ' [W) 

If we now rotate atom 1 about its z-axis by angle 7r, 
then use a third cavity to perform the same filtering op- 
eration on atoms 1 and 2, we will obtain, with probability 
9[-p2 + ( 1 _-p)2]2 the pure GHZ-like state 

I^)i234 = -j= (\egee) + \gegg)) 1234 (17) 

This would, however, require a significantly more compli- 
cated experimental setup. Instead of this final filtering 
stage, we can use an iterative procedure to produce max- 
imally entangled states of arbitrary purity. 

We convert (|16p to a mixture of bipartite entangled 
states by measuring atoms 1 and 2 in the {|±)j} (j = 1, 2) 
basis, where |±)j = (l/v / 2)(|e)j± \g)j), which condition- 
ally produces the outputs 



probability of purification, in Fig.O Obviously, the prob- 
ability of success tends rapidly to zero as the number of 
iterations performed increases to infinity. Although our 
success probability is less than that of Bennett et al., 0] , 
it matches that of Pan et al. [9j . 

Mixed states of the form V'\$+)($+\ + 
(1-7") |*-)(*-| and P'ltf+Xtf+I + (1 - V) |*-)(*-| 
can be purified by means of local conversion into 
the form of Eq. (fl"3| (operated through single-qubit 
rotations) and the application of the scheme depicted 
above. 

In the following two Subsections, we provide two exam- 
ples of applications of the purification protocol we have 
described so far. The objects of our interest are Werner 
states and MEMS. 



Purification of a Werner state 



Let us consider the Werner state 

p w = V\$+)($+\ + Q(l-\$+)($i 



(20) 



where Q = (1 — V)/3. Such states are (for general V) 
mixed, and entangled for V > 4, as witnessed by the 
negative partial transposition criterion. 

In approaching the purification of Eq. ([20]) , some com- 
plications arise due to the fact that cross terms between 
atom pairs in the states |<fr+) and |$_) would not be an- 
nihilated. Indeed, after the first round of purification, 
the output state (conditioned on qubits 1 and 2 being 
projected onto the same state) reads 



P±.34 



7> 2 |$ ± ) 34 ($±| + (1 - V) |*±)34<*±| 



r 2 + {i-ry 



(18) 



The state ' 34 is obtained when atoms 1 and 2 are pro- 
jected onto the same state, and /?_ 34 otherwise. It is 
possible to interconvert p + and p_ by performing a z- 
rotation of angle ir on one of the atoms. This is unnec- 
essary, however, as |3>_) behaves similarly to |<E>+) un- 
der the transformations described above: When applied, 
the protocol will annihilate terms and preserve 

|$_)|$ + ). Likewise, |\&_) behaves analogously to \SSf + ). 

Therefore, if V > |, we can see that we have increased 
the relative proportion of the state |4>±), and hence its 



purity. Correspondingly, Ejj has increased from E 



2V-\toE 



(i) 



(0) 



2V-1 



By iterating this procedure, 



it is possible to produce a state of arbitrary purity. By 
using q replicas of p , the purification protocol produces 
the final state 



v\Q + ) 34 ($ + \ + 0.-v) q \*+) 34 (* + \ 

V q + (1-V) q 



(19) 



with a success probability (P q + (1 - V) q )/2 q . 

The evolution of the input state towards a maximally 
entangled pure state is tracked, in terms of purity and 



p^=b[(v 2 + q")\^ 
+ 4g 2 (|* + )(* 4 



($+| +2^Q|$_)($_ 

+ I*-X*-I)1 



(21) 



with B = V 2 + 2VQ + 5Q 2 . 

This state does not exhibit any great improvement, in 
terms of purity, with respect to pw- Indeed, repeating 





FIG. 5: (a): Purity of the evolved state in the purification 
protocol against the initial purity V and iteration number q. 
(b): Probability of purification against V and q. 
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the procedure using states (|20[) and (|21j) as inputs actu- 
ally causes a decrease in negativity, below that of pw- 

In order to conduct further purification, we prepare 
two atom pairs in state (|2ip , and subject each atom to a 



J rotation around its x-axis. This will effect the trans- 
formations |$ + ) to |$ + ), |$_) to |* + ), |* + ) to |$_) 
and I 1 !'-) to I 1 !'-). Using these two modified pairs as the 
inputs for our protocol, we now obtain the output 



(2) _ [V 4 + 2V 2 Q 2 + 5g 4 ] Ij+Kg+j + 4 [P 2 Q 2 + Q 4 ] (jg-Kg-j + + 8PQ 3 |tt-)(tt_| 

~ V 4 + WV 2 Q 2 + 8V Q 3 + 13Q 4 

I 



which is purified with respect to pw (we can see that it is 
pure to first order in Q). As before, an arbitrarily pure 
state can be obtained by iteration, but the cumulative 
probability of success decreases rapidly with the number 
of steps. 



B. Purification of MEMS 

We now approach MEMS as the second non-trivial 
example of the applicability of our purification proto- 
col. MEMS are states possessing the largest achiev- 
able mixedness for a given degree of entanglement E// > 
[H| ■ A parameterization of MEMS is critically depen- 
dent on the chosen measures of entanglement and pu- 
rity. If Efj is taken to quantify entanglement and the 
linearized entropy Si = (4/3) (1 — Trp 2 ) is used as the 
measure of the purity of a state, MEMS are a single- 
parameter family of states given by 



/ i+yi+3g 2 

6 







V 



2-^/l+3g 2 
3 






Vi+3g 5 



(23) 



with < o < 1 [19j|. MEMS have been theoreti- 
cally characterized and schemes for their generation in 
cavity-QED and circuit-QED have been proposed [33j ] - 
Based on a "Procrustean" method, they have also been 
experimentally produced in an all-optical setup [34j . 
In order to apply the purification protocol, we rewrite 
Pmems. 12 ® Pmems.34 m the Bell basis, let it interact with 
the cavity fields and then exploit the transformations in 
Eqs. (fl~5)) and the analogous transformations valid for the 
remaining tensorial product of two input Bell states. Af- 
ter the projection of qubits 1 and 2 onto the a x eigenbasis, 
we are left with the reduced state of qubits 3 and 4 



o {1) 

rmems, 34 



3^/l+3 g 2 











3^/l+3g 2 




_L 

2 



(24) 




0.4 0.6 0.8 



Lin. Entropy 



FIG. 6: Comparison, in terms of linearized entropy Si and 



negativity E^s, between p„ 



■,34 and p^Ls 



, obtained after 



the purification protocol. The solid (dashed) line shows the 
input (purified) state. The numbered points on each curve 
identify the position of the states associated with eleven values 
of 5 = 0,0.1,..,!. 



This state is purified with respect to the original 
MEMS, and has, for a set value of g, a larger negativity. 
The quantitative comparison between the purified and in- 
put state is shown in Fig. [51 where the states correspond- 
ing to eleven discrete values of g e [0, 1] (equally spaced 
with steps of 0.1) are plotted in the negativity-linear en- 
tropy plane. We recall that, by definition, Si = for pure 
states and Si = 1 for a totally mixed state. Obviously, in 
agreement with the interpretation of MEMS as boundary 
states to physically achievable entangled mixed states, 
the purified states describe a trajectory in the E^f — Si 
plane below the extremal curve traced by proems, 34- 

Evidently, for a set value of g, the corresponding 
Pmems 34 nas a smaller Si , meaning a larger purity, and 
a larger value of negativity, than p m ems- 



C. Simplified purification protocol 

Producing high photon-number Fock states, upon 
which the protocol discussed so far is based, is rather 
difficult. Here, in order to reduce the complexity of 
implementing the scheme, we briefly discuss a simpli- 
fied version in which the two cavities are initially pre- 
pared in single-photon states. The first pair of atoms 
interacts with the cavities for a time t chosen so that 
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Xt = (mi + i)7r, with mi an integer. Analogously, 
the interaction time for the second pair, t' , is such that 
>/2Xt' = (m 2 + i)7r. The evolved states involve cavity 
field modes with up to three photons. 

By projecting the cavity modes onto the single-photon 
state, we realize effective transformations that read, in 
the computational atomic basis 

\eg)i3 -> costfi cos$ 2 |eg)i3 - i(~l) m2 sini?i|#e)i3, 
\ge)i3 -> sini9 2 |eg)i3, 

(25) 

where tfi = ny/2 (mi + §) and i9 2 = ^ (m 2 + |). All 
terms in the input states involving |ee)i3 and |<?g)i3 are 
annihilated by the field-state measurements that, again, 
act as filters. By executing this procedure on two atomic 
pairs, each prepared in state (|13p . we obtain (if atoms 1 
and 2 are both projected onto |+)) the output 

~(i) = V 2 N^){<f'\ + (l-V) 2 N^){^\ 
V 2 N^ + (1 - Vf JV* 

where 

|$'> = == [sin 2 -&2\gg) - sin 2 i?i|ee) + (cosi?i cos?? 2 ) 2 \gg) 
- i(-l) 1 ™ 2 cosi9icost9 2 sini?i {{eg) + \ge))], 

I*') = sintfisintf 2 (|e 5 ) + | 5e )) 

+ cosi9isin(2^ 2 )| 5ff )] 

and i\T$, are normalization factors. By choosing mi 
and m 2 such that cos$i s» cosi9 2 « 0, the states (j2"7) 
reduce to Bell states, and we recover an approximation 
of the ideal n-photon output, (jTHJ). This can be achieved 
with rescaled interaction times of only a few Rabi cy- 
cles, which is advantageous from an experimental point 
of view. For instance, setting m 2 = mi+1 = 3, we obtain 
|($_|$')| 2 = 1-9.8 x 1CT 6 . 

So far, we have overlooked the difficulty of projecting 
cavity modes onto the single-photon state. In our par- 
ticular situation, this is reduced to the simpler problem 
of distinguishing a single-photon state from zero-, two- 
and three-photon states. In fact, the latter of these can 
be effectively removed from the mixture by selecting m 2 
such that sin v2$ 2 ~ (compatible with sin # 2 — 1 for 
m 2 = 7). In order to accomplish the measurement, we 
pass an atom in its ground state through the cavity. The 
interaction time, t" , is selected so that sin (Xt"y/2) = 
and cos (Xt") ss 0. If the cavity was initially in a super- 
position of zero-, one- and two-photon states, the inter- 
action will produce the transformation 

(co|0) + ci|l)+c 2 |2))| 5 )^co|0)| 5 )-ici|0)| e )+c 2 |2)| ff ) 

If, on observing the atom, we find it to be in its excited 
state, we can conclude that the cavity must have been in 
a single-photon state. 



IV. PRACTICAL CONSIDERATIONS 

In this Section, we discuss potential experimental set- 
tings for the implementation of our entanglement con- 
centration and purification schemes. As already men- 
tioned, we will explicitly address microwave cavity- and 
circuit- QED set-ups, which are two of the most promising 
scenarios for quantum information processing and com- 
munication, as stated in the Los Alamos roadmap [35]. 

Microwave cavity-QED uses a mature technology 
based on the interaction between neutral Rydberg atoms 
of large dipole moments and high-finesse microwave cav- 
ities 0, LL3|, |21||. P aradigmaticallv. we follow the proposal 
of Kuhr et al. [l|| , of superconducting cavities sustaining 
a field mode of frequency W//27T = 51 GHz, length L ~ 27 
mm and energy damping time Tu e i4 ~ 130 ms. This cor- 
responds to a cavity finesse F = TTCTfi e id/L ~ 5 x 10 9 , 
where c is the speed of light. As previously mentioned, 
the coherence time of the overall system depends also 
on the damping time of the excited state of the flying 
qubits (each encoded in circular states with principle 
quantum numbers 50 and 51 of a Rydberg atom), which 
is comparable to Tfi e i^. The atomic transition frequency 
(luq = 51.09 GHz) can be tuned in and out of resonance 
with the cavity field by means of a d.c. Stark shift in- 
duced by a static electric field parallel to the cavity axis 
(this field also helps preserve the circularity of the atomic 
states). This technique can be used to control the inter- 
action time between each flying atom and the correspond- 
ing cavity field. The interaction time is also dependent on 
the atomic velocity across the resonator, which is usually 
in the range of hundreds of m.s -1 , with an experimental 
accuracy well- within 5—10%. The effect of this uncertain 
atomic velocity on the concentration protocol is discussed 
in Section[TTl The large atomic dipole moment (~ lOeao 
where e is the electron charge and ao is the Bohr radius) 
allows for the achievement of a strong coupling regime 
characterized by A = 2ir x 50 kHz. This value guar- 
antees the possibility of performing 100 Rabi floppings 
(i.e. Xt ~ 1007t) within r co h e [la ], putting the require- 
ment for a long upper bound on the rescaled interaction 
time (as stated for the concentration protocol) within 
the present state of the art. Atomic state detection, nec- 
essary for both the concentration and the purification 
scheme, can be accurately performed with state-selective 
field-ionization channeltrons (each having a detection- 
error between 10 and 13%). With the new generation 
of high-finesse cavities addressed above, the atomic in- 
jection is performed through large apertures, allowing 
for the screening of the flying qubits from coherence- 
destroying stray fields at the surface of the cavity mirrors. 

In a circuit-QED setup, flying atoms are replaced 
by static superconducting qubits embodied by SQUIDs 
working in the charge regime at the degeneracy point (to 
remove, to first order in the single-Cooper pair charge 
2e, the detrimental effect of low-frequency noise induced 
by background impurities) [13, [23j]. The qubit is inte- 
grated, via conventional optical lithography, in a full- 
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wave, on-chip, coplanar waveguide cavity with resonant 
frequency Uu e id — 5.7 GHz (the microstrip resonator), 
located at a voltage antinode of the sustained field mode 
and capacitively coupled to it. The stripline is a quasi- 
unidimcnsional structure with a very small transversal 
dimension that reduces the effective volume of the cavity 
field and enhances the coupling rate with the qubit. This, 
together with the effective dipole moment of the SQUID 
qubit (~ 2 x 10 4 eao) gives rise to the ratio X/ujfi e id ~ 2%. 
The energy damping time of the stripline is conserva- 
tively assumed to be ~ 1/xs, giving Xt^^ = 100, which 
would, in principle, allow for a large number of coher- 
ent Rabi floppings within the cavity lifetime (experimen- 
tal evidence puts the qubit damping rate in the range 
of 2 /lis), in complete analogy with the cavity-QED set- 
ting described above. A detailed derivation of the qubit- 
stripline coupling Hamiltonian and the resulting coupling 
strength can be found in Paternostro et al. [231 ] . 

The transition energy of the superconducting qubit can 
be adjusted through an external magnetic flux that mod- 
ulates the Josephson energy of the SQUID [12] in such a 
way that the qubit can be easily put in the strong res- 
onant or dispersive regime with the field. This tuning 
ability is the basis of the experimentally demonstrated 
non-demolition measurement of the qubit state, through 
spectroscopic resolution of the field's frequency-pulling 
effect (23|. Therefore, the projection of a qubit required 
by our schemes can be reliably implemented by using an 
ancillary cavity field mode, differing in frequency or po- 
larization with respect to the mode singled out for the re- 
alization of our protocols, along the lines depicted in [23[. 
On the other hand, the same technique can be used in 
order to perform a photon-number-resolving cavity field 
measurement: The qubit spectrum exhibits a peak that 
is shifted, with respect to the empty cavity situation, by 
a quantity depending on the population of the field mode. 

The integration of more than a single SQUID qubit 
in the stripline, in a way that avoids (inductive) cross- 
talk between the qubits, is achievable, the main difficulty 
being the necessity of separating each qubit 's gate volt- 
age [22J to properly set each working point. As opposed 
to cavity-QED, however, the entanglement between su- 
perconducting qubits in circuit-QED has yet to be exper- 
imentally demonstrated, despite promising steps having 
been performed along these lines [361 ]. 

For the sake of completeness, we here mention that a 



scheme for the preparation of MEMS has been suggested, 
both in cavity- and circuit-QED, in Ref. 33]. We refer 
you to these works for an extensive account of the details 
necessary for this procedure. 



V. CONCLUSION 

We have presented two schemes for entanglement con- 
centration and purification for qubit systems. Our pro- 
tocols are explicitly designed in settings involving qubits 
interacting with ancillary cavity fields, which can be im- 
plemented in cavity- and circuit-QED set-ups. Our en- 
tanglement concentration scheme conserves ESPP, there- 
fore being optimal, and achieves the maximum theoreti- 
cal probability of producing a maximally entangled state 
of two qubits. The state purification protocol can be iter- 
ated so as to return a state of arbitrary purity, although 
the success probability decreases exponentially as perfect 
purity is approached. 

The required resources are modest. We need cavities 
prepared in zero- and single-photon states, which have 
been experimentally demonstrated, and the ability to 
perform single-atom rotations and measurements, which 
is possible with high accuracy and large detection efficien- 
cies. In addition, we require cavity-field measurements, 
which can be achieved in an indirect way: We have ad- 
dressed strategies for this step in both of the proposed 
experimental scenarios. 

The mechanisms that constitute the main sources of 
imperfection and error have been quantitatively assessed 
to provide a comprehensive analysis of the protocols at 
hand. We hope that the evidence for state-of-the-art im- 
plementability of our schemes and the astonishing im- 
provements in the control of microwave-based quantum 
technology will pave the way to the realization of en- 
tanglement concentration and purification in cavity- and 
circuit-QED systems. 
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